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PsJ ' Abstract 

> : 

PsJ ■ A simple model, based on the analytical formula for the Compton scattering, 

. is proposed to describe the realistic photon-energy spectra for the Photon Collider 

, at TESLA. Parameters of the model are obtained from the full simulation of the 

I beam by V.Telnov, which includes nonlinear corrections and contributions of higher 

' order processes. Photon energy distribution and polarization, in the high energy 

g_) ' part of the spectra, are well reproduced. Our model can be used for a Monte Carlo 

pi • simulation of gamma-gamma events at various energies and for direct cross-section 

5^ ' calculations. 

> 

u • 1 Introduction 

Photon Collider has been proposed as a natural extension of the e'^e~ linear collider 
project TESLA High-energy photons can be obtained using Compton backscattering 
of the laser light off the high-energy electrons [21 El IH IHl E] • The physics potential of 
the Photon Collider is very rich and complementary to the physics program of e~^e~ and 
hadron colliders It is the ideal place to study the mechanism of the electroweak 
symmetry breaking (EWSB) and properties of the Higgs boson. Precision measurements 
at the Photon Collider may open "new windows" to the physics beyond the Standard 
Model. However, the precise measurements are only possible if the energy spectrum of 
colliding photons is well understood. A detailed simulation of the 77 luminosity spectra 
at the Photon Collider at TESLA has become available recently Pi EI- In this paper a 
simple model based on these results is proposed. 



2 Luminosity spectra 



If the laser beam density is small and the primary electron beams are sufficiently wide, so 
that effects related to the photon scattering angle can be neglected, the energy spectrum 
of the photon beams colliding in the Photon Collider could be calculated directly from 
the Compton scattering cross section (see appendix [X]). However, these assumptions 
will not be fulfilled at the Photon Collider as we need both very powerful lasers and 
strongly focused electron beams to get high luminosity. To find the energy spectra of 
colliding photons, with realistic assumptions about the laser system and the electron 
beams, detailed simulation programs has been prepared [HIE]- 

New samples, with high statistics of the simulated 77 events, have been generated 
recently by V.Telnov 7J. They are based on the beam and laser system layout and 
parameters, as proposed for the Photon Collider at TESLA p]. Electron beam energies 
of 100, 250 and 400 GeV have been considered. Basic parameters used in the simulation 
are listed in Table More details of the simulation can be found in jU E] • 

Shown in Fig.^are the distribution of the colliding photon energy ratio to the primary 
electron beam energy, y = E^/Eg, and the distribution of the 77 center-of-mass energy 
distribution, W^^. Results obtained from the simulation of luminosity spectrum [Zj, for 
electron beam energy of 250 GeV, are compared with the distributions expected from the 
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Table 1: Parameters of the Photon Collider based on TESLA. Listed for different electron 
beam energies Ee are: laser wave length Al, laser photon energy Eq and resulting x param- 
eter values (see section 01 and appendix [X]); horizontal, vertical and longitudinal electron 
bunch sizes a^, ay and a^; average repetition rate, normalized beam emittances •ye^/y, 
/^-functions P^/y , distance between the conversion and interaction point and estimated 
e~e~ luminosity. 
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Figure 1: Energy distribution for photons (left plot) and the 77 center-of-mass energy 
distribution (right plot) from full simulation of luminosity spectrum by V.Telnov [2j (solid 
line), compared to expectations for the simple Compton scattering (dashed line). For 
better comparison of shape, Compton spectra is scaled to the same height of the high 
energy peak. 

simple Compton scattering (lowest order QED), as given by formula (fT^ (see appendix 
IX]) . Realistic simulation indicates that a large fraction of colliding photons will have small 
energies. Also the maximum of the high-energy peak is shifted towards lower energies. 

Shown in Fig. |21 (left plot) is the two-dimensional energy distribution for the colliding 
photons, as obtained from the simulation [Zj, for electron beam energy of 250 GeV. Also 
included in Fig. |21 (right plot) is the energy correlation between two photons calculated as 
the ratio of the two-dimensional energy distribution to the product of two one-dimensional 
energy spectra. Energies of colliding photons are clearly correlated. Majority of collisions 
involve photons with similar energies (large values of the ratio along the diagonal). Col- 
lisions involving one low-energy and one high-energy photon are suppressed (the ratio 
less than 1 in the left-upper and right-lower corner of the plot). This demonstrates that 
the correlation between the angle of Compton backscattering and the photon energy is 
important and has to be taken into account. 



3 Model 

Dedicated studies have been performed in order to understand the differences between 
the photon energy spectra obtained from the collider simulation by V.Telnov [2] and the 
spectra expected for the simple Compton scattering ((T3j). To describe simulation results 
following effects are taken into account: nonlinear effects due to high density of the laser 
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Figure 2: Two-dimensional energy distribution for two colliding photons, obtained from 
the full simulation by V.Telnov [2j (left plot) and the ratio of this distribution to the 
simple product of two one- dimensional energy spectra (right plot). 

beam, correlation between photon energy and the scattering angle, electron rescattering 
and scattering involving two initial state laser photons. 



3.1 Nonlinear effects 

For very high density of the laser beam nonlinear QED effects become important. The 
field of the electromagnetic wave can significantly influence the motion of an electron. 
The effect can be described as an effective increase in the electron mass ml —>■ mg(l + ^^), 
where is the parameter describing the nonlinear effects, proportional to the photon 
density in the laser beam For the Compton scattering nonlinear effects result in an 
effective rescaling of the parameter x, describing the photon energy spectra: 
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where is the electron beam energy and Eq the energy of the laser photon. As a result, 
the energy distribution for the backscattered photons is shifted towards lower energies. 



3.2 Angular correlations 

Electron beams collide with focused laser beams at the distance h ~ 2mm from the 
interaction point. The angular spread of scattered photons is very small (due to very high 
Lorenz boost), but becomes important because of the very small beam spot size. Due to 
the larger scattering angle, the effective vertical size of the photon beam increases at low 
photon energies. As a result, interactions involving low energy photons are suppressed. 
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The effect of angular correlations can be described by the following modification of the 
energy spectrum 

= = /,(,.i).A^exp(-^(i-*-l)). (2) 

where the parameter p relates the vertical beam size ay and the distance b between 
the conversion and interaction points, p ~ {rrieb) / (Egay) , Af is normalization factor and 
fc{y,i) is the Compton spectrum, as described by eq. in appendix 1X1 

3.3 Electron rescattering 

With high density of the laser beam, one can "convert" most of the electrons into the 
high energy photons. However, after the first scattering electrons still have large ener- 
gies. Scattering of laser photons on these secondary electrons results in an additional 
contribution to the photon-energy spectrum. Due to the lower electron energy, secondary 
scattering takes place at lower x value, x' = y'x, where y' is the energy fraction of the 
secondary electrons. Energy distribution for photons scattered off secondary electrons 
can be calculated by integrating over y': 

fiy.i) = ^']dy' wc{y') fc{l-y\S:)-fc{^,y'x) (3) 



where M' is the normalization factor and wc{y') is the weighting function which takes 
into account the dependence of the total Compton scattering cross section on the electron 
energy (SJ • In the high energy limit one gets 

, X \og(x' + 1) 

Photons from scattering on secondary electrons have much "softer" energy spectrum, as 
compared to the Compton scattering on primary electrons. 

3.4 Scattering of two laser photons 

For high density of laser beam it is also possible that electron scatters on two photons 
instead of one. The detailed calculation of the energy distribution for photons produced 
in such scattering is presented in We have found that the distribution obtained from 
the simulation |7j can be well approximated by a simple formula for the scattering on one 
photon with double energy (i.e. with x — > 2x) corrected by an additional factor which 
suppresses the high energy peak: 

^ "^^^ f2{y,x) = /(y,2x).M f--2i-lV , (5) 



N2 dy 



where 5 is the parameter describing suppression of the high energy part of the spectrum 
and ^2 is the normalization constant. 
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4 Parametrization 



4.1 Main assumptions 

The main aim of the presented study was to parametrize the high energy part of the 
luminosity spectra of the Photon Colhder in a simple analytical form. It was assumed 
that the high energy part of the 77 luminosity spectrum can be described as a simple 
product: 

'Tt^~t~ = ftotiyi,x) ftotiy2,x) (6) 

N dyidy2 

where ftotiHi i) is the energy spectrum for the photon. The spectrum can be parametrized 
as a sum of three components described in the previous section: 

ftot{y,x) = nf{y,x) + n f'{y,x) + n2 f2{y,x) (7) 

where n, n' and n2 are parameters describing contributions of different processes to the 
spectrum. All together the model has 10 free parameters which can be adjusted to describe 
the results of simulation by Telnov [7]. Only 4 of these parameters describe the shape of 
the contributing components: 

• two parameters, ^0 and ^1, describing dependence on the electron beam energy: 

e = eO + i?e-6 (8) 

• parameter describing the angular correlations 

• parameter 6 added in description of two-photon scattering 

Remaining 6 parameters, ai . . . ag, are needed to describe the normalization of the con- 
tributing processes: 

• the dependence of the normalization of the Compton scattering contribution [n) on 
the electron beam energy is, in the considered energy range, approximately given 
by: 



n 



M = ai-Ee + as (9) 



• normalization of the electron rescattering {n') and of the two photon scattering (71,2) 
are related to the Compton scattering contribution by the formula: 

n' M' = nAf ■ {as- Ee + a^) (10) 
nsA/'s = nA/"- as (1 - exp(-a6 ^^)) (11) 

Normalization factors A/", A/"' and A/2 are included in the parametrization to simplify 
numerical calculations of the spectra. 

Whereas distributions /, /' and /2 are normalized to unity, normalization of ftot is not 
fixed. It is evaluated from the requirement that the high energy part of the luminosity 
spectrum is given by the formula 
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Figure 3: Comparison of the photon energy distribution obtained from full simulation 
of luminosity spectrum by Telnov [2j , with the fitted contributions of different processes 
considered in the described model, as indicated in the plot. 

4.2 Fit results 

The formula (|7j) was compared with the photon energy spectra obtained from simulation 
by V. Telnov [7 . To minimize effects of energy correlations a cut on the energy of the 
second photon was imposed. For electron beam energy of 100, 250 and 400 GeV the cut 
was 40, 150 and 260 GeV, respectively. Parameters of the model were fitted to the photon 
spectra, for y > 0.1, simultaneously at all energies. 

Result of the fit to the photon energy distribution at Ee=250 GeV is shown in Fig. El 
Fitted contributions of different processes are also indicated. The model describes the 
spectra very well down to -E^ ~ 0.1 -Eg. Three processes considered in the model contribute 
to different parts of the spectra. By summing these contributions most of details of 
the distribution can be well reproduced. Very good description of the photon energy 
distribution is obtained for all considered energies, as shown in Fig. |3 

Normalization of the fitted parametrization, as well as of the contributions of different 
processes, are shown in Fig.Elas a function of the electron beam energy Eg. Normalization 
of the parametrization changes from about 0.8 at 50 GeV to about 0.55 at 500 GeV. This 
means that the two photon spectrum obtained from the product of the two distributions, 
as given by eq. describes between 65% and 30% of events expected from the spectra 
simulation [7j. 35% to 70% of the total 77 luminosity expected from simulation is due to 
events with one or two low energy photons not described by our parametrization. 

In Fig. ini the comparison of the average photon polarization resulting from the fitted 
parametrization with the distribution obtained from the simulation of luminosity spectra 
is shown. To describe the photon polarization two additional assumptions were made in 
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Figure 4: Comparison of the photon energy distribution from the fitted parametrization 
with the distribution obtained from full simulation of luminosity spectra for three 
electron beam energies, as indicated in the plot. Imposed cut on the energy of the second 
photon is 40, 150 and 260 GeV respectively. 
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Figure 5: Normalization of the CompAZ parametrization of the photon energy distribu- 
tion, relative to the distribution obtained from full simulation of luminosity spectrum [7] , 
as a function of the electron beam energy. Also shown are normalizations of separate 
processes considered in the model. 
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Figure 6: Comparison of the photon polarization resulting from the fitted parametrization 
with the distribution obtained from full simulation of luminosity spectra [Zj, for three 
electron beam energies, as indicated in the plot. Imposed cut on the energy of the second 
photon is 40, 150 and 260 GeV respectively. 

the model: 

• scattering involving two photons results in very high photon polarization. It is taken 
from the Compton formula (fTTj) (with x 2x) for scattered photon energies above 
the threshold for one photon scattering, and fixed at the threshold value for lower 
energies; 

• electrons undergoing secondary scattering are unpolarized. 

Both assumptions have no strong physical motivation,^ however they were found to give 
the best description at the simulation level. It has to be stressed that the model was 
not fitted to the photon polarization distribution and no additional parameters were 
introduced to describe it. Very good agreement between the parametrization and the 
average photon polarization obtained from the simulation, is observed for y > 0.3. 

4.3 CompAZ 

The routine implementing the described spectra parametrization is called CompAZ. It 
can be used to calculate the photon energy spectrum for different electron beam energies 
and the average photon polarization for a given photon energy. Separate contributions 
from three considered processes ()2I3I5|) can also be calculated. Additional routines were 

^As it was pointed out by V.Telnov, significant polarization is expected for secondary electrons |12) . 
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Figure 7: Comparison of the center-of-mass energy distribution obtained with the Com- 
pAZ parametrization with the distribution obtained from full simulation of luminosity 
spectra [Zj, for three electron beam energies, as indicated in the plot. 

prepared for convenient event generation from the parametrized spectrum. All routines 
can be downloaded from web P^ . 

In Fig.[7|the 77 center-of-mass energy distribution obtained from CompAZ is compared 
with the distribution obtained from the simulation of luminosity spectra [Zj , for different 
electron beam energies. No cuts on photon energies were imposed. Proper description of 
the spectra is obtained for W^^ >~ 0.3 Wmax, where Wmax = 2 E^ax is the maximum 
center-of-mass energy available for two photons produced in the Compton scattering. Also 
the average product of the photon polarizations, related to the ratio of 77 collisions with 
the total angular momentum = and \Jz\ = 2, is properly described for large W^^ 
(see Fig. |H1). Center-of-mass energy distribution for two colliding photons with = 
is shown in Fig. IHl The parametrization describes very well the high energy part of the 
spectra, most relevant for many physics studies. 

In the calculation of the 77 luminosity spectrum as a product of two energy distri- 
butions (jni) possible energy correlations between two beams are neglected. In Fig. ^1 
the two-dimensional energy distribution obtained from CompAZ, and the ratio of this 
distribution to the two-photon spectrum obtained from the simulation by V.Telnov [Zj 
are shown. In the high energy part of the spectrum, when both photons have high ener- 
gies, y > 0.5 Umax, this ratio is close to 1. This shows that CompAZ properly describes 
this part of the spectrum and no additional corrections for beam energy correlations are 
needed. Only for W^^ <~ 0.3 Wmax energy correlations become important. CompAZ 
overestimates the number of collisions with one high energy {y > 0.5 ymax) and one low 
energy {y < 0.3 ymax) photon (the ratio greater than 1), and underestimates the number 
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Figure 8: Comparison of the average product of photon polarizations from the CompAZ 
parametrization with the distribution obtained from full simulation of luminosity spectra 
[7j, for three electron beam energies, as indicated in the plot. 
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Figure 9: Center-of-mass energy distribution for two colliding photons with = 0. 
Results obtained with the CompAZ parametrization are compared with the distribution 
obtained from full simulation of luminosity spectra ^ , for three electron beam energies, 
as indicated in the plot. 
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Figure 10: Two-dimensional energy distribution for two colliding photons, obtained from 
the CompAZ parametrization (left plot) and the ratio of this distribution to the one 
obtained from full simulation of luminosity spectrum (right plot). 

of collisions involving two low energy photons (the ratio smaller than 1). 
4.4 Applications 

As already mentioned in the previous section, dedicated routines are available for fast sim- 
ulation of 77 scattering events with CompAZ. They have been recently used in the sim- 
ulation of W^W~ and ZZ pair-production at the Photon Collider, for different electron 
beam energies Distribution of the 77 center-of-mass energy, W^^ for 77 W^W~ 
events generated with PYTHIA, for electron beam energy of 250 GeV, is shown in Fig.lTTl 
Generated events were reweighted for photon polarizations. Sample of events generated 
using CompAZ is compared with the sample generated with the luminosity spectrum from 
simulation [HI EI- Very good agreement is observed. The advantage of CompAZ is that it 
can be easily used for any beam energy^, giving reasonable description of the energy and 
polarization. 

CompAZ parametrization can also be used for calculating expected event distributions 
without the time consuming event generation. Numerical integration is few orders of 
magnitude faster than the full simulation and can be used to extrapolate results of full 
simulation to other beam energies. Results from the recent study of the heavy Higgs 
boson production at the Photon Collider for Higgs mass of 180 GeV {h — > ly+VT") 
and 300 GeV {h ZZ) are shown in Fig. ^| Expected invariant mass distributions 
obtained from the full simulation, based on PYTHIA J3j and fast detector simulation 
program SIMDET 16^ , are compared with results obtained by the numerical convolution 

^Parametrization can be used for 50 < i?e < 500 GeV. 
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Figure 11: Distribution of the 77 center-of-mass energy for 77 W^W~ events generated 
with PYTHIA, for electron beam energy of 250 GeV. Sample of events generated using 
the CompAZ parametrization is compared with the sample generated with the luminosity 
spectrum from full simulation. Generated events were reweighted for photon polarizations. 

of the cross section formula for vector boson production in 77 scattering with the CompAZ 
photon energy spectrum and parametrization of the detector resolution. The agreement 
between both approaches is very good. With numerical integration based on CompAZ, it 
was possible to calculate the detector level effects expected from the interference between 
direct W^W~ production and h —> W~^W~ decay. To estimate the effect with full event 
and detector simulation very large statistic of events would be required. 



5 Summary 

Luminosity spectrum obtained from the detailed beam simulation is the best tool for ac- 
curate simulation of 77 interactions. Recently the new version of CIRCE code became 
available ^3 QB] , which contains the Photon Collider luminosity spectra based on simu- 
lation by V.Telnov [HI EI- The CIRCE program [IB] gives detailed description of the lumi- 
nosity spectra in the whole energy range taking properly into account all non-factorizing 
contributions (energy and polarization correlations). The package includes routines for 
convenient event generation. However, only three selected electron beam energies have 
been considered so far {Ee= 100, 250 and 400 GeV). Moreover, parameters of the Photon 
Collider assumed in the simulation are only known with accuracy up to 10-20%, as many 
details of the project are still not fixed. Therefore, other models resulting in the similar 
(or better) accuracy are also applicable for detailed studies. 

We propose the model describing the photon energy spectra of the Photon Collider 
at TESLA in a simple analytical form, based on the formula for the Compton scattering. 
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Figure 12: Distribution of the reconstructed invariant mass for 77 — > W'^W~ events with 
a SM Higgs-boson mass of 180 GcV and an electron-beam energy of 152.5 GeV (left plot) 
and for 77 — > ZZ events, with a SM Higgs-boson mass of 300 GeV and an electron- 
beam energy of 250 GcV (right plot). Results from the simulation based on PYTHIA 
and on the SIMDET detector simulation (histogram) are compared with the distribution 
obtained by the numerical convolution of the cross-section formula with the CompAZ 
photon energy spectra and parametrization of the detector resolution (solid line). The 
distribution expected without the Higgs contribution is also shown (dashed line). 

Parameters of the model are obtained from the comparison with the full beam simula- 
tion by V.Telnov, which includes nonlinear corrections and contributions of higher order 
processes. Photon energy distribution and polarization, in the high energy part of the 
spectra, are well reproduced in a wide range of electron beam energies. Model can be 
used for Monte Carlo simulation of gamma-gamma events. Parametrization is also very 
useful for numerical cross-section calculations. 



Acknowledgments 

Special thanks are due to V.Telnov for providing the results of the Photon Collider lumi- 
nosity spectra simulation and to M.Krawczyk for many productive discussions and valu- 
able comments. I would also hke to thank V.Telnov, I.Ginzburg and T.Ohl for comments 
and critical remarks to this paper. 



14 



A Energy spectrum for Compton scattering 



The energy spectrum of the photons resulting from the Compton backscattering of laser 
light off the high energy electron beam depends on the electron beam and laser polariza- 
tions, Pe and Pl, and on the dimension-less parameter x: 
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(12) 



where Eg is the electron beam energy, Eq the energy of the laser photon and rrie is the 
electron mass. The maximum energy of the scattered photon is 
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and the energy spectrum is given by ^ 
1 dN 
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where r = y/{x{l — y)), y is the fraction of the electron energy transfered to the photon 
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and Afc is the normalization factor given by 
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The degree of polarization of the photons scattered with energy fraction y is given by 
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(17) 



The energy spectrum and polarization of the scattered photons, for a; = 4.5 and 85% 
polarization of the electron beam (proposed parameters of the TESLA Photon Collider 
for electron beam energy i?e=250 GeV), and various helicities of laser beam are shown in 
Fig.d 
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Figure 13: Energy distribution (left plot) and polarization (right plot) for photons from 
Compton back-scattering, for different laser beam polarizations Pl, as indicated in the 
plot. X = 4.5 corresponds to laser wave length of 1.06 fim and primary electron beam 
energy of 250 GeV. Electron beam polarization Pe=85%. 
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